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I. INTRODUCTION 

Most systems in nature cannot be regarded as purely closed systems but open ones, interacting with an environment 
(e.g. a thermal bath). The evolution of these systems must then be non unitary, i.e., interactions with the environment 
must lead to dissipation as well stochastic effects, which is the way the environment backreacts on the system. One 
common way for describing this non unitary evolution is by means of stochastic Langcvin-likc equations of motion. 
These non-deterministic equations of motion are used in many systems of interest, such as simulating Brownian motion 
in (classical and quantum) statistical mechanics and in other areas of physical interest [l[ . 

A classical example of a system whose dynamics is modeled by a Langevin equation of motion is the one that 
describes the Brownian motion of a classical particle of coordinate q, unitary mass and subjected to a potential V{q) 
(as usual dots mean derivative with respect to time and y[(7(i)] = dV/dq), 

qit)+vm + v'm]^m, (i-i) 

where is a Markovian (local) dissipation term and ^{t) is a stochastic term with white noise and Gaussian properties, 
satisfying (throughout this work we consider the Boltzmann constant /cs = 1) 

mm)^2Trjs{t-t'). (1.2) 

Approaches with Langevin equations such as Eq. (jl.ip and its generalizations, are used in different contexts, e.g. 
in classical statistical mechanics to study problems with dissipation and noise, to determine how order parameters 
equilibrate and in the studies such as dynamic scaling and dynamic critical phenomena [2, 3]. 

Though extensively used, equations of the form of Eq. with noise properties as given by Eq. (jl.2p . can only be 

considered phenomenologically. This is because it implicitly assumes that the environment interacts instantaneously 
with the system. This is a physically unacceptable situation that violates causality, since the environment bath has 
no memory time. Its is worth mentioning that a similar situation happens with the description of the dynamics of 
a conserved order parameter by a Cahn-Hilliard equation which lacks causality f5|. In the case of the Cahn- 
Hilliard equation this happens because, since it is a diffusion-reaction type of equation, it should be characterized 
by microscopic scattering events. In real systems scattering events proceed through finite time intervals, which, 
consequently, must lead to finite memory effects. In order to fix this bad behavior of the Cahn-Hilliard equation, 
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memory effects must be taken into account, as explicitly shown recently in Ref. [bI]. In a microscopic description of the 
effects of the environment degrees of freedom on some select variables taken as the system, the same is expected to 
happen. Dissipation and stochastic (noise) terms are expected to originate from scattering events, thus giving origin 
to finite interaction times, that reflect in the system's equation of motion as nonlocal (i.e. non-Markovian) terms 
with memory effects. The simplest archetype of this is the description of the system-environment as being modeled 
by linearly coupled harmonic oscillators [7] (for a general review, see e.g. Ref. 8]), which also become to be known as 
Caldeira-Leggett type of models [9|]. The derived equation of motion for the system variable, when the bath degrees 
of freedom are integrated out, leads to a generalized Langevin equation (GLE) of the form 

m + f dt'D{t - t') q{t') + V' [q{t)] = m , (1-3) 

where <o is some initial time, D{t — t') is a dissipation kernel and the noise term £^{t) is still Gaussian with zero mean 
but colored, i.e., with two-point correlation, according to the Fluctuation-Dissipation theorem for classical systems, 
and, thus £,{t) satisfies 

mm) = TD{t^t'), (1.4) 

where averages are assumed to be taken with respect to a bath of free oscillators with equilibrium distribution p^' 
at some temperature T. Similar derivations in the context of field theory models (see, for instance Refs. [13: [HI) also 
show the emergence of generalized equations of motion of the form of Eq. (|1.3p . with a more complicated structure 
depending on the form of the coupling of the system (e.g. some field we are interested in the dynamics) with the 
environment (or bath fields, made of the remaining fields other than the one representing the system field). 

Despite the fundamental differences between Eqs. (jl.ip and (|1.3p . we expect that when the time scale of the memory 
kernel D{t~t') is much smaller than any other time scales of the system, the local approximation (jl.ip . with dissipation 
coefficient defined by [l^l 

r/ = / dt'D{t - t') , (1.5) 

can still be a good approximation for the system's dynamics and the effects of finite memory be negligible. There is 
also an immense saving of effort as well as much more transparent understanding of the physics from a local equation 
as opposed to a nonlocal one, since the former can generally be analyzed with much less numerical treatment than 
the latter, thus it is a very important question to know when, and how accurately, a generally nonlocal equation 
can be approximated by a local form. Likewise, when the memory constrains cannot be ignored, such as when the 
typical microscopic time scales are large in comparison to the other time scales characterizing the dynamics, we must 
then be able to have appropriate tools to tackle the non-Markovian equations of motion. This is possible with some 
restrict forms of kernels, as we are going to see below, which, nevertheless, can represent physically relevant systems 
of interest, thus, being well motivated. 

In this work our objective is to gauge the applicability of an equation of motion of the form of Eq. (jl.ip when 
compared to the non-Markovian form, when considering some of the most common forms for the dissipation kernel 
D{t,t'). These forms for the dissipation kernel include for example the one that describes an Ornstein-Uhlenbeck 
(OU) process [l3| and the exponential damped harmonic (EDH) kernel [l^ (see also Ref. [l^ for a recent review 
on the different colored noise terms used in the literature). We study both the cases of additive and multiplicative 
noises, including system dependent dissipation terms, according to the Fluctuation-Dissipation theorem. A detailed 
numerical analysis is made when the various parameters characterizing the thermal bath, e.g. the bath relaxation (or 
damping) parameter, frequency and temperature of the bath are varied. 

The remaining of this work is organized as follows. In Sec. II we define the prescription to transform the non- 
Markovian equations in a system of Markovian time differential equations. We study specifically the OU and EDH 
kernels. In Sec. Ill we present the numerical results for the non-Markovian equations for different model parameters 
and compare the results to those coming from their Markovian approximations. The results are obtained for both 
the cases of additive and multiplicative noises. How the dynamics depends on the various parameters characterizing 
the thermal bath is studied in details. In Sec. IV we present our conclusions, discuss the various results we have 
obtained, and we give a possible relation and the relevance of our results for the study of nonequilibrium dynamics 
in field theory models. 
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II. THE NON-MARKOVIAN LANGEVIN-LIKE EQUATION 

Here we study a GLE describing the interaction of a system, denoted by a variable (j) (which can be e.g. the 
coordinate of a particle) in interaction with a thermal bath, where the noise has the properties such as in Eq. (|1.4p . 
The GLE studied here is then of the generic form 

4>{t) + r{t) f dt'cf>^\t')D{t - 1') m + v'{4>) = <^"(t) m , (2.1) 



where n = 0, 1, with n = giving the standard GLE of additive form, such as Eq. (|1.3p . while for n = 1 gives a 
multiplicative GLE. The multiplicative noise and system dependent dissipation form are motivated from field theory 
calculations [lOj. ill] and this is why we have also included this special case here in view of future applications in that 
case. The potential in Eq. (j2.ip is considered to be one with quadratic and quartic terms, given by 

n0) = ^0^ + ^0^ (2.2) 

where and A are parameters depending on the details of the system under study. Here, we can associate m with 
the system's frequency and A with the degree of nonlinearity of the system's potential (or in the context of field 
theories, with the strength of the system's self-interactions). 

Nonlinear GLEs of the form of Eq. (|2.ip are notably difficult to solve. Analytical methods can only be used when 
the equation can be approximated or put in a linear form, such as in the additive noise case and when the quartic 
term in the system's potential can be neglected. This is because, in the additive noise and variable independent 
dissipation case, such as in Eq. (II. 3p . the equation is in the form of a convolution, so can be solved through Laplace 
transform for instance Otherwise, in the more general cases, we must resort to numerical methods. This is the 
approach we follow in this paper in order to analyze the dynamics obtained from Eq. (j2.ip . Though there are some 
specific numerical methods using e.g. Fourier transform that may apply for equations with non-Markovian kernels 
of generic form [l^, we still would like to be able to solve equations such as Eq. (|2.ip through standard methods, 
which are less numerically expensive than other alternatives. This is the case, for example, when using Runge-Kutta 
methods. Recently, the authors [17| have demonstrated the reliability of using a fourth-order Runge-Kutta method 
when solving GLE of the OU and EDH forms. The way this can be done stems from the fact that non-Markovian 
equations with kernels of those forms can be replaced by a system of completely local first-order differential equations, 
which has been described in details in p7j . 

As already mentioned, in this work we concentrate our study in equations such as Eq. (j2.ip with non-Markovian 
kernels of either the OU type 



Dou{t-t') = rjje-'<<^'-''\ (2.3) 



or with the EDH type [Ti 



Dnit-t') = ^e-'^(*-*')|t{cos[ni(t-i')] 

+ ^sin[f]i(t-i')]| , (2.4) 

where in the equations above, rj sets the magnitude of the dissipation, 7 sets the relaxation time for the bath kernels, 
T = 1/7, and r^o gives the oscillation time scale in the case of the EDH kernel. In Eq. (|2.4p . 17^ = 17q — 7^, and so, in 
the EDH case the values of 7 and fig are restricted so to have ilf > 0. 

It can be easily shown that the OU and EDH noises can be generated by the stationary part of the solution of the 
following differential equations, respectively, 

iou{t) ^ -l[^ou{t) - V^m] , (2.5) 

^H{t) + 2jiH{t) + nl^Hit) = nl^W^at) , (2.6) 



4 



where C in Eqs. (12. 5p and (|2.6[) is a white Gaussian noise satisfying 



(CWC(i')) = Sit-t'). (2.7) 



Taking to = 0, we can define a new variable u'(t) by [T3, [T^ 



t 

w;(t) = - / dt'(j,'^{t')D{t~t')^{t') . (2.8) 





This, together with Eqs. (|2.5p and (|2.6[) . leads to the following system of local first-order differential equations 
representing the GLE Eq. ()2.1|1 : For the OU case, 

= y , 

y = -l/'(0) + 0"w;oc/ + <^"eoc/ , 
Co;/ = 

while for the EDH case we obtain 



£,ou-V^C , (2.9) 



= 2/, 

2/ = -y'(0)+0"u;^f + (/."Cff , 

u = + 1)^(0) 0"j/_27i?^^(O)0"y, 



i = + , (2.10) 

where the variable u{t) in Eq. (|2.10p is defined as 



u{t)= I dt' 







dDnit- t') 



dt' 



nt')'-^ , (2.11) 



and in Eqs. ((^ and ([^11})) . we also have that l?oc/(0) = ?7 7, £'//(0) = 77f^g/(27) and 1)^(0) = 0, which follow from 
Eqs. (|2.3p and (|2.4p . The additive noise case is when n = is taken in Eqs. (|2.9p and (|2.10p . while ri = 1 is for the 
multiplicative noise case. These two cases, for both types of kernels, are next numerically studied below. 

III. CONTRASTING THE NON-MARKOVIAN NUMERICAL RESULTS WITH THE MARKOVIAN 

APPROXIMATION ONES 

Let us now consider our numerical results for the Markovian and non-Markovian dynamics for the system. In the 
Markovian approximation all memory effects are neglected and the non-Markovian dissipation term in Eq. (j2.ip is 
replaced by a local dissipation term with magnitude as given by Eq. (|1.5p . i.e., we write Eq. (|2.ip in the form 

0(t) + V 0'"(i) Ht) + ^'(0) = nt) m • (3.1) 

As we have mentioned before, in general we expect the local form the GLE to be a valid approximation when the 
relaxation time scale for the thermal bath, r = I/7, is much smaller than the characteristic time scale for the system, 
e.g., T <^ 4>/4> (this is equivalent to the quasi-adiabatic condition set in Ref. [l^ in the field theory case for the validity 
of the local Markovian approximation). When this condition is met in a sufficiently large time interval At = t — ipj 
thus At/r ^ 1 (which is equivalent as taking to ~* —00) and the time nonlocal term in Eq. (j2.ip can be written as 
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f dt'D{t-t')(j)''{t')^{t') 

J to 

f -t')^r] (j)^"{t) ^{t) , (3.2) 

J to — > — QO 

where in the last step we have used the definition Eq. (|1.5|) . The result p.2p then leads to the local dissipation 
term in Eq. (|3.ip . Of course, under the conditions set above, at sufficiently short times we expect the memory 
effects to influence the dynamics in some significant way, but these memory effects should quickly become negligible 
at long times, t 3> I/7. In any event, we thus expect that after some long time period the memory effects can 
become sufficiently damped such that the Markovian approximation could represent well the overall dynamics of 
the system. After all, we expect that both dynamics, the non-Markovian and the Markovian ones to both have the 
same asymptotic state. But we still face with a natural and important question to answer: For a given set of model 
parameters representing the system and the thermal bath to which it is coupled to, for how long can we expect 
the memory effects due to the non-Markovian terms to be important and when can they be neglected and then the 
dynamics be well represented by the local Eq. (|3.ip ? This is because, even though both dynamics are expected to 
approach each other asymptotically, the time this happens could be so long that the memory effects could lead to 
important physical effects and the local Markovian dynamics would just not be appropriate to be used. Since the 
representation of the dynamics in a local form as given by Eq. (j3.ip represents a considerable simplification, for both 
a numerical point of view, or for analytical analysis (when it is possible), when compared, e.g., to the full nonlocal, 
integro-differential stochastic Eq. (|2.ip . this then becomes an important question to be accessed for most practical 
studies that make use of stochastic equations of motion. It is also important to investigate how the dynamics is 
affected by varying not only 7, but the other parameters characterizing the thermal bath, such as the temperature T 
and frequency fio, whose effects on the dynamics may not so direct as the ones obtained by just varying 7. Below we 
try to answer all these questions, performing our study, numerically, in the context of the additive and multiplicative 
noise cases, with either the OU or EDH kernel terms, defined in the previous section. 

Next we show the results of our systematic simulations of the system of differential first-order equations, Eqs. (|2.9p 
and (|2.10p . for the non-Markovian GLE with OU and EDH kernels, respectively. The results are compared to those 
obtained through the local approximation given by Eq. (|3.ip . All our simulations were performed with 300 000 
realizations over the noise and we have integrated all differential equations using a standard fourth-order Runge- 
Kutta method with a time stepsize varying between St = 0.01 and St = 0.001, which were found to be more than 
enough for both numerical stability and also for enough numerical precision (as determined in Ref. jl7j . these values 
already assure an overall numerical error of always smaller than about one percent, which suffices for our comparison 
purposes set here). In all our simulations we have also used the initial conditions 0(0) = 1 and 0(0) = 0. The 
time in all our evolutions is in units of the (inverse of the) frequency for the system (which is equivalent to consider 
m = 1 throughout). Comparisons between the Markovian and non-Markovian dynamics are made varying the relevant 
parameters of the bath for the two cases of kernels considered, while keeping the system parameters fixed. 

A. The additive noise case 

Let us now turn to our numerical results. The relevant bath parameters are the dissipation magnitude 77, the 
temperature T (that are common to both Markovian and non-Markovian dynamics), the bath damping parameter 7, 
and the bath frequency Qq (in the EDH kernel case). Since the dissipation magnitude is common for both types of 
dynamics, in the following we keep 77 fixed, at the value rj ~ \ throughout (which can be checked to correspond to lead 
to an underdamped dynamics in the local equations of motion) and we vary the remaining bath parameters. This will 
allow us to better understand the importance of the memory effects alone for the dynamics. So, we consider the various 
dynamics as 7, JIq and T are changed. We include the study of how the dynamics changes with the temperature 
because this is useful to determine how this internal property of the thermal bath influences the dynamics when 
comparing to both the Markovian and non-Markovian cases. 

We start our analysis by first considering variations in the parameter 7. Representative values for 7 are then chosen 
and all the remaining parameters are initially kept fixed. Note that since 7 acts as damping the effects of the nonlocal 
kernels, the larger is 7 the better must be the local approximation for the full nonlocal dynamics. Then we keep 7 
fixed at the largest value used in our analysis and then consider variations in the other parameters. This will allow us 
to determine the consequent importance of the remaining parameters and whether a change of those parameters can 
discriminate the types of dynamics studied, for example discriminate additive and multiplicative stochastic dynamics. 

Let us first consider the GLE with additive noise and OU kernel. This is considered in Fig. [1] where we plot side 
by side our results for the dynamics of the ensemble averaged macroscopic system variable 0, (0) = </3(i), where the 
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Figure 1: OU case with additive noise and its IVlarkovian approximation: time evolution for 'p{t). (a) 7 = 0.5, (b) 7 — 1.0 and (c) 
7 = 5.0. The other parameters are taken as Qo ~ 1.0, rj — 1.0, T = 1.0, m = 1.0 and A = 1.0. 



average is over the noise realizations. This is obtained from the Markovian and non-Markovian equations, Eqs. (|3.ip 
and (|2.9p . respectively, with n = 0. 

In Fig. [2] we plot side by side our results for (p{t) for the Markovian and non-Markovian regimes for the EDH case, 
again by considering the additive noise (n — 0). 

The effect of changing 7 seen in both Figs. [T]and[2]is clear and well within the expected: The larger is the relaxation 
time (1/7) for the nonlocal kernels, the larger are the memory effects, resulting in a strong difference with respect to 
the local approximation, seen most notably at short times. As also expected, at some sufficient long time, that we 
here see to depend on how large 7 is, the two dynamics, Markovian and non-Markovian approximate one of the other. 
This is also seen if we had plotted the correlation {4>'^{t)) for both OU and EDH cases. The difference between the 
Markovian and non-Markovian dynamics can also be better estimated by defining the quantities, 



non— Markovian Markovian 
non— Markovian /Markovian ■ 

(3.3) 

The results for the differences A(p and are shown in Figs. [3]and[4l for the OU and EDH cases, respectively. 
The results from the plots shown in Figs. [3] and [3] are useful to determine within which time scale the Markovian and 
non-Markovian dynamics become sufficiently close (within to some given precision) . The results for these time scales 
for the different simulations we have performed for the Markovian dynamics and for the non-Markovian dynamics 
with the two types of memory kernels explored here will be given below in Table [H 

An immediate conclusion we can realize from an inspection of the results shown in Figs. [T]-|3]is that the time 
scale that it takes for the two dynamics to begin to be equivalent is much larger than the time scale for the kernel 
relaxation itself, I/7, and also larger than the typical system's time scale, which is typically given by the inverse of 
the system's frequency (1/m). This will remain true for the multiplicative noise case with either the OU or the EDH 
memory kernels. 
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Figure 2: EDH case with additive noise and the IVlarkovian approximation: time evolution for <^(t). (a) for 7 — 0.1, (b) for 7 = 0.3 
and (c) for 7 = 0.5. The other parameters are the same as in Fig. 1. 
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Figure 3: The differences A(j> and A^^ in the OU additive case. All other parameter kept fixed as before. 

B. The multiplicative noise case 

Let us now verify the results concerning the Markovian and non-Markovian stochastic equations when multiphcative 
noise and system dependent dissipation are concerned. We then here explore the case n = 1 in Eqs. (|2.9p and (|2.10p . 
for the OU and EDH cases, respectively. 

In Fig. Owe plot side by side our results for {(j>) for the Markovian and non-Markovian regimes in the OU case with 
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Figure 4: The differences A</> and A^^ in the difference in the EDH addit ive case. 
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Figure 5: OU case with multiplicative noise and its Markovian approximation: time evolution for ip{t). (a) for 7 — 0.5, (b) for 
7 = 1.0 and (c) for 7 = 5.0. The other parameters are taken as S7o ~ 1-0, rj — 1.0, T = 1.0 and A — 1.0. 

n = 1, while in Fig. [6] are the results for the EDH case. 

The results for the differences and A(/)^, for the n = 1 multiplicative noise case, are shown in Figs. [7]and[8l for 
the OU and EDH cases, respectively. 

The results shown in Figs. [5] -[5] again indicate that, as expected and similar to the additive noise n = case, that 
the two dynamics, local and nonlocal, become closer to each other the larger is the kernel damping parameter. They 
also seem to indicate that in the multiplicative noise case the two dynamics are closer to each other for the same 
parameters used in the additive noise case. This can be more quantitatively estimated through the differences ()3.3|) . 
The results for the time scales when the two dynamics begin to be sufficiently close to each other are also presented 
in Table n 

We note from the results shown in Table [J that the system variable (p in the case of the GLE with OU additive noise 
tends to approximate the corresponding Markovian dynamics faster than in the case of the multiplicative noise case. 
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Figure 6: EDH case with multiplicative noise and the Markovian approximation; time evolution for ip{t). (a) for 7 = 0.1, (b) for 
7 = 0.3 and (c) for 7 = 0.5. All other parameters taken as in the previous figures. 
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Figure 7: The differences A0 and A</)^ for the OU multiplicative case. All other parameter kept fixed as before. 



The opposite seems to happen to the equal time correlation function {(jy^), where in the case of OU multiplicative 
noise it is faster than in the case of additive noise. The behavior of A(p in the additive noise cases indicates that 

the averaged system variable {(j)) is approached faster to the Markovian dynamics than in the multiplicative noise 
cases. The behavior for the dynamics of (0^), when looking at the behavior of A^^, is analogous, except in the EDH 
multiplicative noise case, where it is slower than in the additive noise case. 

In addition to the above results obtained for the ensemble average system variable (j), it is also useful to determine 
how the memory effects influence the thermalization time for the system when put in contact with the thermal bath 



10 




0.2 
0.1 
j- 0^ 
-0.1 
-0.2 



Y=0.1 
Y=0.3 
y=0.5 









5 10 


15 

t 


20 25 


30 


5 


10 15 20 

t 


25 30 






Figure 8 


: The differences A<j) an 


d A(?!>^ 


for the EDH multiplicative case. 




7 


^'t'OV add 


^"^OU add 


A0OU mult 


mult 


7 


A't'EDH add 


add 


A(^EDH mult 


^<^EDH mult 


0.5 


73 


28 


35 


70 


0.1 


180 


89 


30 


73 


1.0 


32 


25 


30 


24 


0.3 


67 


33 


28 


71 


5.0 


14 


7 


27 


12 


0.5 


53 


26 


23 


67 



Table I: The approximate time scale, in units of 1/m, for the non-Markovian dynamics to approach the Markovian one, within 



a precision of 10 for the differences defined in (|3.3|l . 



at some temperature T. For this, let us define an effective time dependent temperature for the system according to 
the equipartition of kinetic energy. 



Teff(i) = {^\t)) . (3.4) 

The results for T^s for the OU and EDH cases for the additive and multiplicative noise situations are shown in Fig. 
[9l where it is also plotted the Markovian, local approximation case, for comparison. 

From the plots shown in Fig. [9] we again see the overall behavior seen in the previous plots. The larger is the 
relaxation time for the memory kernels (i.e., the smaller is the damping parameter 7 of the memory kernels), the 
larger the dynamics takes to approach the Markovian approximated one. We also see how memory effects reflects in 
the thermalization of the system. Larger bath relaxation times lead to a longer time for the system to thermalize. 
Typically, for comparable bath relaxation time scales, in the additive case the system tends to thermalize faster than 
in the multiplicative case. In Table [TTl we give the approximate time (in units of 1/m) for thermalization for all 
the cases studied above and where this behavior for the time for thermalization in each of the types of dynamics 
analyzed can be verified. Note that in all cases, with additive or multiplicative noise, the Markovian dynamics tend in 
general to underestimate the time scale for thermalization, when compared to the non-Markovian dynamics, except 
for large 7 in the OU additive noise case, where the thermalization in the non-Markovian dynamics tends to be better 
approximated by the Markovian one. As far the different dynamics are concerned, we observe from both the plots in 
Fig. [9] as from the results shown in Table [TTl for the thermalization times, that the additive noise case always tend to 
have a smaller thermalization time than the multiplicative noise case. This is observed for both the Markovian and 
non-Markovian cases. 



C. Markovian and non-Markovian dynamics in terms of the temperature of the thermal bath 

We now explore how the temperature of the thermal bath will influence the results shown in the previous section. 
We again concentrate on the differences between the Markovian and non-Markovian dynamics and the thermalization 
time. For this study we consider the cases with the highest values of 7 considered in the previous subsection, which 
gives the best comparison between the dynamics. Keeping the highest values of 7 allows us to determine whether the 
comparison between the two dynamics improves or worsens as the temperature is changed. 

The results for the differences Acj) and A^^ for the GLE with OU kernel are shown in Fig. [TUl while for the GLE 
with EDH kernel are shown in Fig. [TlJ The time scale for which the full non-Markovian dynamics approaches the 
respective Markovian dynamic approximation, in each of the cases studied here, are tabulated in Table IHII 

From the results shown in Table IHII we can see that higher temperatures for the thermal bath give only minor 
improvements in the OU additive noise case as regarding the approach of the non-Markovian dynamics to the ap- 
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Figure 9; The effective temperature Te// for the (a) OU n — case, (b) OU n = 1 case, (c) EDH n = case and (d) EDH n = 1 
case. 
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Table II: The approximate time for thermalization, in units of 1/m, for the Markovian and non-Markovian dynamics, determined 
when l|3.4|l approaches the temperature of the thermal bath. 



proximated Markovian one. But in the multiplicative (OU) case the changes are much stronger, with (</)) and (0^) 
approaching much faster to the Markovian dynamics the larger is the temperature. In the EDH case, for the additive 
noise, ((/>) tends to approach slower to the Markovian approximation and in the multiplicative noise case the approach 
is much faster, compared to the additive one. But the correlation (0^) changes much differently. In the EDH additive 
noise case the Markovian approximation worsens as the temperature is increased, while in the multiplicative case the 
Markovian approximation improves, but only slightly for high temperatures. Similar behavior to this will also be 
seen below for the thermalization times for each of the dynamics. We can also observe from the results shown in the 
plots of Figs. [TOl and fTTl that for the multiplicative noise case, for both non-Markovian dynamics studied, the memory 
effects are much stronger at lower temperatures than at high temperatures and these effects last for a much longer 
time than in the additive noise cases. 

We now study how the thermalization time for each of the dynamics studied here will change when the temperature 
is changed. In Fig. [12] we show the plots for the various cases of dynamics studied here. The thermalization times 
for each case are tabulated in Table IIVI 

From the results shown in the plots of Fig. [12] and in Table IIV| we can clearly see that while the cases with OU 
noise the thermalization times for the Markovian and non-Markovian dynamics are similar (using the largest values 
of 7 considered in Table [T| , the non-Markovian dynamics with EDH noise, for both the additive and multiplicative 
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Figure 10: The results for and A(^^ for the OU n — case, (a) and (b), respectively, and for the OU n — 1 case, (c) and (d), 
respectively. 



T 


A<^OU add 


^4>OV add 


^<?^OU mult 


^4>OV mult 


At/^EDH add 


^4>EDH add 


A(/>EDH mult 


^4>E'DH mult 


0.1 


17 


10 


140 


206 


68 


24 


126 


220 


1.0 


14 


7 


27 


12 


53 


26 


23 


67 


10.0 


13 


6 


6 


7 


42 


49 


13 


61 



Table IIL The approximate time scale, in units of 1/m, for the non-Markovian dynamics to approach the Markovian one, 
within a precision of 10~^ for the differences defined in (|3.3p for different temperatures for the thermal bath. The parameter 7 
considered is the largest one considered in Table |T] in the OU (7 = 5.0) and EDH cases (7 — 0.5). 

cases, produces thermalization times considerably larger than those of the Markovian cases. As in the case of Ac/y^ for 
the EDH additive noise case seen in Table IIIIl we again see here the anomalous behavior of the thermalization time 
increasing as the temperature is increased in the EDH additive noise case. This is opposite to the behavior seen in 
the OU and EDH noise cases with multiplicative noise. Note also that in the case closer to the Markovian dynamics, 
the OU additive noise case, there is almost no relevant difference in thermalization times as the temperature is varied 
up to 2 orders of magnitude. For other values of temperature of the thermal bath we tested, we found that these 
conclusions do not change. 

D. The non-Markovian dynamics in terms of the frequency of the thermal bath 

We now investigate how the frequency of the thermal bath affects the dynamics of the system. Note that in this 
case only the GLE with EDH kernel is affected (the Markovian and OU dynamics are independent of the frequency 
of the thermal bath). As before, we can better interpret the results by analyzing the differences (|3.3p and the graph 
of the effective temperature. From them we obtain the time scale for the non-Markovian dynamics to approach the 
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Figure 11: The results for and A^^ for the EDH n — Q case, (a) and (b), and for the EDH n — 1 case, (c) and (d). 
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Table IV: The approximate time for thermalization, in units of 1/m, for the Markovian and non-Markovian dynamics, deter- 
mined when (|3.4|l approaches the temperature of the thermal bath. 



Markovian one and we can better estimate the differences in thermalization times (if any) in both cases. The results 
for this case, when the frequency of thermal bath is varied, are shown in Fig. [T3]for the differences (|3.3p . while the 
behavior of the effective temperature, when the frequency of the bath is changed, is shown in Fig. [TH In all cases 
the frequency Oq is in units of the frequency of the system, m, and the time is in units of 1/m. As in the previous 
analysis for the EDH case, we consider 7 = 0.5, and all other parameters kept fixed at the values as given before, 
except fio- The time scales for the non-Markovian and Markovian dynamics to approach to each other are given in 
Table |Vl while the time for thermalization for each case is given in Table IVII 

From the results seen in Fig. [13] and Table IVl we can note that the smaller is the frequency of the thermal bath the 
larger it takes for both {(f)) and (^^) in the non-Markovian case to approach the Markovian approximation. As the 
frequency of the bath is increased beyond the frequency of the system, the time scale for the non-Markovian dynamics 
to approach the Markovian one tends to decrease, but it also rapidly reaches a point where larger frequencies for the 
bath would not make the Markovian approximation to improve too much. 

The results are seen to be much different again when the thermalization of the system is studied. From the 
thermalization plots seen in Fig. [THand the thermalization times given in Table IVTl we see that smaller frequencies 
for the bath, compared to the system's frequency, lead to much higher thermalization times for the non-Markovian 
approximation compared to the Markovian one. The thermalization time for the non-Markovian dynamics is seen 
to decrease with an increase of fio, except the multiplicative noise case, where above a certain point ^ "^j the 
thermalization time starts to increase again. Though the importance of the memory kernel is expected to be less 
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Figure 12: The results for the effective temperature as a function of time, when the temperature of the thermal bath is changed, for 
the cases of (a) Markovian additive noise, (b) IVlarl<ovian multiplicative noise, (c) OU additive noise, (d) OU multiplicative noise, (e) 
EDH additive noise and (f) EDH multiplicative noise. 



important the larger is flo, which makes it fast oscillate, the behavior in the multipHcative noise is probably far less 
simple, because the higher are the nonlinearities in that case, which are brought by the system dependent dissipation 
and noise. 
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Figure 13: The results for and A(^^, when the frequency of the thermal bath (in units of m) is changed, for the EDH n = 
case, plots (a) and (b), respectively, and for the EDH n = 1 case, plots (c) and (d), respectively. 




t t 

Figure 14: The results for the effective temperature as a function of time, when the frequency of the thermal bath (in units of m) 
is changed, for the cases of (a) EDH additive noise and (b) EDH multiplicative noise. 

IV. CONCLUSIONS 

In this work we have analyzed in details the differences between the dynamics of a system when treating it in terms 
of its full non-Markovian equation of motion and when expressing it in terms of its Markovian, or local approximated 
form. Having set the appropriate description for the non-Markovian equations, we have then studied the applicability 
of the local approximation for these equations. Wc here have concentrated in two forms for the non-Markovian memory 
kernel: The OU and EDH cases and we have analyzed the cases of additive and multiplicative noises in both cases. 
We have seen that in general, for most of the parameters in both cases, the local approximation is far from being to 
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Table V: The approximate time scale, in units of 1/m, for the non-Markovian dynamics to approach the Markovian one, within 
a precision of 10^"' for the differences defined in p.3p for different frequencies of the thermal bath. The parameter 7 considered 
is the largest one considered in Table |T] in the EDH cases (7 = 0.5). 
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Table VI: The approximate time for thermalization, in units of 1/m, for different frequencies of the thermal bath in the in the 
Markovian approximation and in the EDH additive and multiplicative noise cases. 

represent a good description of the dynamics. Obviously, since these are all dissipative systems, we expect the two 
dynamics, non-Markovian and Markovian, to tend to each other asymptotically. We have then analyzed how long it 
takes for each of the non-Markovian dynamics studied to tend to the Markovian ones. We have analyzed both the time 
for the system variable, (</>), as well for the equal time correlation, {(j)^)- We have also analyzed the thermalization 
time for each of the dynamics by studying the behavior of the correlatfon function {(jP), which, according to the 
equipartition theorem, can be associated to the temperature of equilibration of the system. 

Each of the dynamics were investigated by changing the main parameters characterizing the thermal bath: The 
damping term for the memory kernels, 7, the temperature T and the frequency of the thermal bath. The 
parameters of the system were kept fixed for convenience as well the magnitude of the dissipation, 77, which is a linear 
parameter entering in all the dynamics and that was kept fixed at the point where all the dynamics were initially 
underdamped. By increasing 7, the memory kernels are damped faster and the Markovian approximation tends to be 
better as expected. As 7 is varied, besides the expected behavior of the non-Markovian dynamics to approach faster 
the Markovian one the larger is the 7, it is also observed that for the same value for the memory kernel damping term, 
in general the cases with EDH non-Markovian dynamics tend to approach faster to the Markovian dynamics than the 
OU cases, as far the thermalization times in each of the dynamics are concerned. It is also observed that in general 
the Markovian dynamics tend to overestimate the time for thermalization as compared to the time it takes in the 
non-Markovian cases. The thermalization times in the studied dynamics also tend to be larger in the multiplicative 
noise cases than in the additive ones. 

When analyzing the behavior of each of the dynamics when varying the temperature and frequency of the thermal 
bath, we have fixed the value of 7 and then investigated each of the time scales for non-Markovian dynamics for ((/>) 
and (0^) to approach the corresponding Markovian approximation. By increasing the temperature it is observed that 
the Markovian approximation tends to improve, except in the case of the dynamics of the correlation (0^) in the EDH 
case with additive noise, where the Markovian approximation tends to worsen and in the OU additive noise case, 
where the dynamics is weakly dependent on the temperature of the thermal bath. These results seem to indicate that 
the effective dissipation in the other dynamics is dependent on the temperature of the thermal bath, in particular in 
the multiplicative noise cases, where the effective dissipation seems to increase much faster with the temperature than 
in the additive noise cases. The thermalization times shown in Table HVl are also indicative of this behavior. Finally, in 
the study of how the frequency of the thermal bath affects the dynamics (in the EDH memory kernel cases) , we have 
seen that as the frequency of the thermal bath is increased, the Markovian approximation seems to improve till some 
frequency close to the system's frequency. Above that value there is little improvement. We have also verified that 
for frequencies of the thermal bath much below the system's one, the Markovian approximation worsens considerably. 

A few generic results can also be drawn from the analysis of all cases studied here. In particular, we can note 
from the obtained results that either the local approximation underestimates the effective dissipation seen in the non- 
Markovian dynamics, or overestimates it in most of the regions of parameters. The local (Markovian) approximation 
for the dynamics tends to be better at larger values of the bath damping term 7 and for larger values of the frequency 
Hq and temperature T of the thermal bath (except for the correlation {(j)'^) in the case of EDH memory kernel with 
additive noise). The difference between the dynamics is larger at short times, exactly as expected because of the finite 
memory times for the non-Markovian equations. The different simulations we have performed with different bath 
parameters allowed us to estimate the approximate time scales when the Markovian approximation may become an 
appropriate description of the dynamics. In general, this time scale is much larger in the multiplicative noise cases 
than in the additive noise ones. Also, given the specific differences seen in each of the dynamics when varying e.g. 
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either 7, T or Ho, by looking at these differences when changing the properties of the thermal bath may be a useful 
way to discriminate possible stochastic phenomena in nature and to tell whether they can be dominated by additive 
or multiplicative noises. 

Finally, we should point out the possible connection and relevance of our studies with those related to the dynamics 
of field theory models. It has been shown [l^, [ll[ that non-Markovian kernels of the same form as studied here 
can also appear in the studies of the effective dynamics of an order parameter in field theories and in cosmology 
in general. Since the studies of dissipative processes in those applications are related to time non-local terms in 
the effective evolution equation of the system, for example, for a background scalar field describing the system, or 
an order parameter for a phase transition problem, we expect our results to be of relevance for understanding the 
relevant dynamics in those situations as well. In particular, our studies may be useful to clarify the applicability or 
not of approximating the dynamics in those problems as local ones, as usually it is considered to be the case there. 
The results we have obtained here shows that, in many cases, the local approximation is not a reliable description of 
the true non-Markovian dynamics. The difference can be very large at short times and continue to be for long time 
scales, with memory effects making a strong contribution for the dynamics. This may have strong consequences, for 
example, when studying thermalization and equilibration times in phase transition problems, or in the problem of the 
production of particles and radiation in cosmology. In the dynamics of some systems in contact with a thermal bath, 
as we have seen in the studies performed in this work, the usual local Langevin equation typically underestimates the 
thermalization with respect to the true dynamics, indicating that the use of local approximated forms for the study 
of the dynamics can be unappropriated and even lead to erroneous results as regarding to the system's equilibration 
and thermalization time scales. 
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